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Abstract 



We consider two-grid mixed-finite element schemes for tfie spatial discretiza- 
tion of the incompressible Navier-Stokes equations. A standard mixed-finite 
element method is applied over the coarse grid to approximate the nonlinear 
Navier-Stokes equations while a linear evolutionary problem is solved over 
the fine grid. The previously computed Galerkin approximation to the ve- 
locity is used to linearize the convective term. For the analysis we take into 
account the lack of regularity of the solutions of the Navier-Stokes equations 
at the initial time in the absence of nonlocal compatibility conditions of the 
data. Optimal error bounds are obtained. 



1 Introduction 



In this paper we study two-grid mixed finite-element (MFE) methods for 
the spatial discretization of the incompressible Navier-Stokes equations 

ut- Au+ {u- V)n + \7p = /, (1) 
div(u) = 0, 

in a bounded domain J7 C M"^ (d = 2, 3) with a smooth boundary subject to 
homogeneous Dirichlet boundary conditions u = on 90. In ([T]), n is the 
velocity field, p the pressure, and / a given force field. As in [23], [21], [25] 
we assume that the fluid density and viscosity have been normalized by an 
adequate change of scale in space and time. We approximate the solution u 
and p corresponding to a given initial condition 

u{-,0)=uo. (2) 

Two-grid methods are a well established technique for nonlinear steady prob- 
lems, see [33]) [35] . The main idea in a two- level method involves the dis- 
cretization of the equations over two meshes of different size. A nonlinear 
system over the coarse mesh is solved in the first step of the method. In 
a second step, a linearized equation based on the approximation over the 
coarse mesh is solved on the fine mesh. In [28], |29j several two-level meth- 
ods are considered to approximate the steady Navier-Stokes equations. In 
these papers, depending on the algorithm, the second step is based on the 
solution of a discrete Stokes problem, a linear Oseen problem or one step of 
the Newton method over the fine mesh with the coarse mesh approximation 
as initial guess. 

Several two-level or two-grid schemes have also been considered in the lit- 
erature to approximate the evolutionary nonlinear Navier-Stokes equations 
([TJ-(l2|). Again, two approximations to the velocity (and correspondingly 
two approximations to the pressure), are computed. One of them is defined 
by a discretization of the nonlinear equations over a coarse mesh, uh, and 
another one, Uh, is defined by an appropriate linearization over a fine mesh. 
Different classes of algorithms can be seen as two level methods. In par- 
ticular, although they were originally developed from different ideas, the so 
called nonlinear Galerkin methods, postprocessed and dynamical postpro- 
cessed methods, fall into this category. 

Postprocessed Galerkin methods were first introduced for Fourier spec- 
tral methods in [18], [19] (see also |31j ) and later extended to finite element 
methods in [7] , [6] , [H] , [15] . In all these works the main idea is the following: 
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one first compute the standard Galerkin approximation to the velocity and 
pressure over a coarse mesh {uh,Ph) of size H and then compute the post- 
processed approximation in a finer mesh at selected times in which one wants 
to obtain an improved approximation. More precisely, the postprocessed ap- 
proximation {uh^Ph) computed at a given time t* is an approximation in a 
mesh of size h <^ H to the following (steady) Stokes problem: 

-Au + Vp = f- iuH{t*) - {unit*) • V)uH{t*)) 

div(n) = Oj 

u = 0, 

Here, unit), t G (0, T], is the standard MFE approximation computed in the 
coarse mesh in a time interval (0,T] and t* £ (0,T]. Note that the compu- 
tation of {uH(t),pH{t)), t G (0,T], is completely independent of the compu- 
tation of {uh{t*),ph(t*)) in the fine mesh. The postprocessed approximation 
improves the rate of convergence of the standard Galerkin approximation 
over the coarse mesh in the following sense. If the rate of convergence of 
the Galerkin approximation to the velocity in the L^(r2)'^ (j = 0) or H^{Q)'^ 
(j = 1) norm is 0{H^~^) then the rate of convergence of the postprocessed 
approximation to the velocity is 0{H^^^~'^ \log[H)\)+0{h^'~^). Analogous 
results are obtained for the pressure. For first order mixed finite element 
methods the improvement in the rate of convergence of the velocity is only 
achieved in the H^^QY norm, [6j. Then, if one wants to achieve the opti- 
mal accuracy of the fine level in the H^{^})'^ norm, one can first compute 
the Galerkin approximation on a coarse mesh of size H = and then 

compute the postprocessed approximation over the fine mesh of size h at the 
desired time levels. For example, one should take H = /i^a and H = /i^/s for 
linear and quadratic mixed finite elements, respectively. It can be expected 
that the computational cost of the postprocessed approximation is smaller 
than that of the Galerkin approximation on the same fine mesh, since in the 
former method the time evolution is done on the coarse mesh, and only at 
selected time levels are computations done on the fine mesh. This has been 
confirmed by the numerical experiments in [7] (see also jl6j and jl8j ) 

In [31j a related algorithm, the so-called dynamical postprocessing, is 
introduced for the Fourier case. In this algorithm, the standard Galerkin 
approximation, {uh,Ph), is computed over a coarse mesh in the first level, 
as before. For the second level an approximation to a linear evolutionary 
problem, instead of the steady problem is computed. More precisely, 
the dynamical postprocessing involves the approximation, at each time step. 
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over a mesh of size h <^ H of the problem: 

f^u-Au + Vp = f-{uH-V)uH\ . ^ 

> m 
div(n) = Oj (4) 

u = 0, on do,. 

Note that in the dynamical postprocessing, the computation of {uH{t)-,PH{t)) 
and {uh{t),ph{t)), t G [0, T], is coupled. The rate of convergence of the dy- 
namical postprocessing scheme is proved in [31] to be the same as the rate of 
convergence of the standard postprocessing. In the case of highly oscillatory 
solutions the dynamical algorithm is shown to be more efficient than the 
standard postprocessing in some one dimensional examples. The dynamical 
postprocessing method is also considered in [33], named now as two-level 
method, in the case of mixed finite elements. In [33], the author treats the 
fully discrete case integrating in time with the backward Euler method. A 
similar two-level scheme is also considered and analyzed in [21J where the 
author uses first order mixed finite elements in space, Crank-Nicolson ex- 
trapolation for the time integration over the coarse mesh and the backward 
Euler method for the time integration over the fine mesh. 

The so-called nonlinear Galerkin methods are also two-level methods 
that have been used to compute approximations to ([I])-©. They were first 
introduced for Fourier spectral methods [12], [32], and later extended to 
mixed finite element methods in [5]. In this work the authors obtain the 
rate of convergence of the nonlinear Galerkin method in the case of first or- 
der elements. The rate of convergence is the same one of the postprocessed 
method. The main difference between the nonlinear Galerkin methods and 
the postprocessed or two-grid methods is that in the former the approxima- 
tion on the coarse mesh takes into account the influence of the fine mesh, 
whereas in the latter it is just the standard Galerkin method (i.e., computed 
independently of the fine mesh). 

In this paper we analyze two two-grid algorithms in the context of spatial 
mixed finite element discretizations to approximate the solutions of ([H)-©. 
The two algorithms we consider are very similar to the dynamical postpro- 
cessing method. The difference is the treatment of the nonlinearity in the 
second level. In the dynamical postprocessing method the nonlinear con- 
vective term of the fine level is approximated by the data {uh ■ ^)uh (see 
the right-hand-side of (|4])). In the two algorithms we consider in the present 
paper, the approximation to the velocity of the coarse mesh uh is used to 
linearize the nonlinear convective term of the fine level. In the first algo- 
rithm, the linearized convective term of the fine level is {uh ■ V)uh- In the 
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second algorithm uh is regarded as an initial guess to perform one Newton 
step in the fine level. For the spatial discretization we consider mixed fi- 
nite elements of first, second and third order. More precisely, we consider 
the mini-element and the quadratic and cubic Hood- Taylor elements. The 
analysis for other mixed finite elements of the same order is similar. As in 
[21], [m due to the lack of regularity at t = of the solution of ([I])-© no 
better than 0{H^) error bounds can be expected. For this reason we do not 
analyze higher than cubic finite element discretizations. For the temporal 
discretization we use the backward Euler method or the two-step backward 
differentiation formula. The analysis of the fully discrete methods is similar 
to the one appeared in [15] and it is only briefly indicated in this paper. 

This is not the first time these two algorithms have been considered. The 
first algorithm was introduced in [20], where the authors analyze the semi- 
discrete in space case for first order finite elements. In [2] the authors extend 
this analysis to the fully discrete case and in [1] the second order Hood- Taylor 
finite element is used for the spatial discretization and the two-step backward 
differentiation formula for the time integration. In [27] the second algorithm 
is analyzed for the Fourier spectral case while in [30] the analysis is extended 
to the case of first order mixed finite elements considering the fully discrete 
case coupled with the Crank-Nicolson scheme for the time integration. As 
opposed to the above mentioned works on the same methods, in the present 
paper we take into account the lack of regularity suffered by the solutions of 
the Navier-Stokes equations at the initial time. Then, for the analysis in the 
present paper we do not assume more than second-order spatial derivatives 
bounded in up to initial time t = 0, since demanding further regularity 
requires the data to satisfy nonlocal compatibility conditions unlikely to be 
fulfilled in practical situations [23] , [23] . This is the first time these methods 
are analyzed under realistic regularity assumptions. Also, this is the first 
time the cubic case is considered and the first time the quadratic case is 
considered for the second method. 

There are some other improvements with respect to previous works. In 
[1] the authors get an error bound of order 0{H^ + + (At)^) for the 
fine approximation to the velocity Uh in the H^{D.)'^ norm whenever the 
following inequality is satisfied aiH"^ < (At)^ < a2H^, ai and 02 being 
constants independent of H and At. With the technique of this paper 
an error bound of order 0(| log(/i)|| log(i/)|-ff^ + h'^ + (At)^) for the same 
fully discrete method in the H^{Q.)'^ norm can be obtained for H and At 
independently chosen. With the new error bound obtained in this paper one 
can achieve the rate of convergence of the fine mesh in the H^[Q)'^ norm by 
taking H = h^/'^ instead oi H = h"^/^. This fact improves the efficiency of the 
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method compared with the (same order) standard Galerkin method over the 
fine mesh. Also, the authors of [I] remark that they have observed the same 
rate of convergence for the two-grid method with H = h^^"^ and H = h?^^ in 
the numerical tests they have carried out, which supports the improved rate 
of convergence we obtain in this paper. We want to remark that in all the 
numerical experiments of |27] . |30) and [1] the two-grid algorithms improve 
the efficiency of the standard Galerkin method in the sense that a given error 
can be achieved with less computational cost with the new algorithms than 
with the standard Galerkin method. In |27) a comparison in the Fourier case 
between the standard postprocessing, the dynamical postprocessing and the 
second two-grid algorithm is also included. Although the computational 
cost of the two-grid approximation over the fine mesh is bigger than that of 
the postprocessed approximations, the two-grid algorithm produces smaller 
errors in the case of moderate to high Reynolds numbers. Finally, comparing 
the two algorithms we analyze in this paper we remark that with the second 
algorithm better error bounds are obtained in terms of H. Although this fact 
could make the choice of the second algorithm preferable for computations, 
it turns out in practice to be rather inefficient to solve the linear systems 
accurately. For this reason, some authors suggest solving instead an Oseen 
problem leading then to the first algorithm, see [28j . 

The rest of the paper is as follows. In Section 2 we introduce some 
preliminaries and notations. In Section 3 we carry out the error analysis of 
the first two-grid algorithm in the semi-discrete in space case. In Section 
4 we consider the analysis of the second two-grid algorithm in the semi- 
discrete in space case. Finally, in Section 5 we consider the fully discrete 
case integrating in time with the backward Euler method or the two-step 
backward differentiation formula. 

2 Preliminaries and notations 

We will assume that i7 is a bounded domain in M*^, d = 2,3, not necessarily 
convex and smooth enough. When dealing with linear elements (r = 2 
below) Q may also be a convex polygonal or polyhedral domain. We will 
consider the Hilbert spaces 

F = {n e {L\n)f I div(u) = 0,u- n^^^ = 0}, 
V = {u£ {H^{n) f I div(u) = 0}, 

endowed with the inner product of L^(r2)'^ and Hq{Q,)'^, respectively. For I > 
integer and 1 < q < oo, we consider the standard Sobolev spaces W^''^{Q)'^ 
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of functions with derivatives up to order / in and H\^Y = W^''^{^Y. 

We will denote by || • \\i the norm in H^{^)'^, and || • ||_/ will represent the 
norm of its dual space. We consider also the quotient spaces //'(J7)/M with 
norm \\p\\hi/m. = inf{||p + c||; | c G M}. 

Let us recall the following Sobolev's imbedding inequalities [IJ: For q G 
[1,00), there exists a constant C = C{VL,q) such that 

lis 

(5) 

For g' = 00, di]) holds with | < f ■ 

Let n : L'^{9.Y — > H he the L'^{9.Y projection onto H. We denote by 
A the Stokes operator on 

A:V{A)C1H — >H, A = -HA, V{A) = H'^i^f r\V. 

Applying Leray's projector to ([T|), the equations can be written in the form 

ut + Au + B{u,u) = Ilf in 0, 

where B{u,v) = n(n • V)f for u, v in Hq{Q)'^. 

We shall use the trilinear form •, •) defined by 

b{u,v,w) = {F{u,v),w) \/u,v,w £ -ffQ($7)'^, 

where 

F{u, v) = {u- V)v + ^(V • u)v Vn, v G H^{nf. 

It is straightforward to verify that b enjoys the skew-symmetry property 

b{u,v,w) = —b{u,w,v) \/u,v,w G HQ{il.)'^ . (6) 

Let us observe that B{u,v) = IlF{u,v) for u e V, v G HQ{^y. 
We shall assume that 

< Ml, ||n(t)||2 < M2, 0<t<r, 

and, for k >2 integer, 

L(fc-2)/2j 

sup ||4''''^V||fc_l_2|fc/2| + Y] sup < +00, 
7=0 
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so that, according to Theorems 2.4 and 2.5 in [23j . there exist positive 
constants Mfc and K/^ such that for k > 2 

\Ht)h + \\Mmk-2 + Mmm-^/R < MkT{tf-^'^ (7) 

and for A; > 3 

j <yk~z{s){\\u{s)\\l^\Us{s)\\l_2M\p{s)\\]jk~l^ < KI, 

(8) 

where T(t) = min(t, 1) and (7„ = e '*V"'(s) for some a > 0. Observe 
that, for t < T < oo, we can take r(t) = t and crn{s) = s". For simphcity, 
we wiU take these values of r and cj„. We note that no further than k < 6 
will be needed in the present paper. 

Let Th = {t^, 4>^)i£i^ , /i > 0, be a family of partitions of suitable domains 
il/i, where h is the maximum diameter of the elements r/* G Th and 0^ are 
the mappings of the reference simplex tq onto t^. We restrict ourselves to 
quasi- uniform and regular meshes Th- 

Let r > 2, we consider the finite-element spaces 

Sh,r = {xh(^C{Tlh)\Xh\r^O(t>'i e P''-\to)], 

Sir = {Xh G C{Tlh) I Xh\r^ o 0^^^ G ^''"'(^o), Xh{x) = Vx G aj],,}, 

where P^'~^(to) denotes the space of polynomials of degree at most r — 1 on 
tq. Since we are assuming that the meshes are quasi-uniform, the following 
inverse inequality holds for each Vh G {Sy^^Y (see, e.g., llO,, Theorem 3.2.6]) 

\\vh\\w"^>i{TY - C^'"""'^^^"' V/tllvK'.9'(r)''' (9) 

where 0</<m< 1, 1 < q' < q < oo, and r is an element in the partition 
Th. 

We shall denote by {Xh^riQh,r-i) the so-called Hood-Taylor element 
[H [26], when r > 3, where 

Xh,r = {Slrf , Qh,r^l = Sh,r-^1 □ ^'(0^)/^, r > 3, 

and the so-called mini-element P] when r = 2, where Qh,i = 'S'/i,2nL^(0/j)/R, 
and Xh^2 = {S^ 2)'^ ® ^h- Here, Mh is spanned by the bubble functions br, 
T G Th, defined by br{x) = {d + l)'^'^^Ai(3;) • • • Arf+i(x), if x G r and else- 
where, where Xi{x) , . . . , Xd+i{x) denote the barycentric coordinates of x. 
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For these mixed elements a uniform inf-sup condition is satisfied (see [8]); 
that is, there exists a constant /3 > independent of the mesh grid size h 
such that 

. r (Qh,^ ■ Vh) . r, 

mf sup J — 7^ — 7 > p. (10) 

qheQh,r-l Vh€Xh,r \\'"h\\i\\qh\\L^ 

The approximate velocity belongs to the discretely divergence-free space 

Vh,r = Xh,r n {Xh G HliVih) I (^ft, V • Xh) =0 Vg/, G Qh-r-l]- 

We observe that, for the Hood-Taylor element, Vh^r is not a subspace of V . 
Let II/j : L?{0Y — ^ Vh,r be the discrete Leray's projection defined by 

(n^u, Xh) = {u, Xh) yxh G Vh,r- 

We will use the following well known bounds 

||(/-n;,)u||^. < C/i'-J'||n||,, l<l<r, j = 0,l. (11) 
Assuming that is has a smooth enough boundary, we also have 

j|^-m/2jj(^_ jj^^^ll^ < ^^/+min(m,r-2) ii^ll^^ 1 < / < r, m = 1, 2. 

(12) 

Since {A^^^'^Hhf ,Vh) = (/,^^^^^Wh), for all Vh G V/,,,., it follows that 

p-^/^n,/||o < (13) 

Moreover it holds for / e L^($7)'^, see [2]: 

llV^'n^/llo < Ch'\\fh + \\A-''^Ilf\\o s = l,2. (14) 

Let A denote either A = A oi A = Ah- Notice that both are positive self- 
adjoint operators with compact resolvent in H and Vh, respectively. Let us 
consider then for a G M and t > the operators A°^ and e"*-^, which are 
defined by means of the spectral properties of A (see, e.g., [ll, p. 33], p!7]). 
An easy calculation shows that 

P"e-*-^||o < (ae-^)"t-", a > 0, t > 0, (15) 

where, here and in what follows, ||-||q when applied to an operator denotes 
the associated operator norm. Also, using the change of variables t = s/t, 
it is easy to show that 

{»-"l<V<-.-l|„..,^B(l,i), ,10, 

where B is the Beta function (see, e.g., [3j). 
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3 Semi-discretization in space. The first two-grid 
algorithm. 

In this section we carry out the error analysis of the first two-grid algorithm 
for the Hood- Taylor mixed finite element with r = 3 or 4. At the end of 
the section we include the results that can be obtained for the mini-element 
with a similar but simpler analysis than the one showed along the section. 

The first algorithm we consider is the following. Let us choose h < H 
so that Vn^r C Vh,r- Then, in the first level we compute the standard mixed 
finite-element approximation to ([I])-©. That is, given uh{0) = IIh{uo), 
we compute unit) G Xh^^ and pnit) G QH,r-i, t G (0,T], satisfying, for all 

e XH,r and TpH G QH,r^l 

{uhAh) + {VuH,y(l)H) + Huh, uh^cPh) + {"^Ph, <Ph) = (/, <Ph) (17) 

(V ■UH,iJH)=0. (18) 

In the second level we solve a linearized problem on a finer grid and given 
UhiO) = UhUo we compute Uh{t) G Xh,r and ph{t) G Qh,r-i, t G (0,r], 
satisfying, for all 4>h G X^^r and iljh G Qh,r-i 

{uh, 4>h) + (Vn/i, V0fe) + {uH • Vuh, (t>h) + {^Ph, 4>h) = (/, <Ph) (19) 

{V ■Uh,i^h) = 0. (20) 

To obtain the error bounds for {uh,ph) we will follow the error analysis of jl4^ 
and introduce an auxiliary approximation (see [14[ Section 4.1]). For a u and 
p solution of ([I])-© let us consider the approximations Vh : [0, T] — > X^ ^ 
and gh : [0,T] — > Qh,r-i, respectively, solutions of 

{vh, <Ph) + (Vfh, \/(ph) + [Vgh, M = if, <Ph) - Ku, u, <ph) (21) 
{V ■Vh,i^h) = 0, (22) 

for all (ph £ Xh r and iph G Q/i.r-i; with initial condition Vh{0) = IlhU-o- We 
will also use the following notation: 

Zh = UhU - Vh- (23) 

Next, we state some lemmas that are needed in the proof of the main theo- 
rems. The first one summarizes previous results. 
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Lemma 1 Let {u,p) be the solution of ([H)-©. There exists a positive con- 
stant C such that the discrete velocity defined by \21\) -[2 ^) and the Hood- 
Taylor element approximation to u, uh, satisfy the following hounds for 
j = 0,1, andt e (0,r]; 

\\vH{t)-UH{t)\\j < j^^\\og{H)\H'+'-\ 3<r<4, (24) 

\\uH{t)-u{t)\\, < -^-^H'-^, 2<r<5, (25) 

\\uh{s) -u{s)fjds < CH^^^-^\ 3<r<4. (26) 

Proof The bound (j24|) is proved in Theorems 4.7 and 4.15 in [T3]. The 
case j = in (j26p is proved in [23l Theorem 3.1] and [241 Theorem 3.1]. 
The case j = 1 follows from the case j = by applying and p!T|) . see 
also Corollaries 4.8 and 4.16 in [T3]. Finally, (j26p is proved in Lemmas 5.1 
and 5.2 in p4|. □ 

For the convenience of the reader, we will reproduce here the following two 
Lemmas, the first one from [7] and the second one from 



Lemma 2 For any f € C([0, T]; L^(0)'^), the following estimate holds for 
allte [0,T]: 

|^^e-(t-M.n,/(s)||od5 < C|log(/i)| max^||/(t)||o. 

Lemma 3 Let {u,p) be the solution o/ ([T])~([2]). Then, there exists a positive 
constant C such that the error Zh = H/iU — Vh in ( t^3|) satisfies the following 
hound: 

\\A[^'^'^"z,,\\,<-^^h^+'-=, j = 0,1,2, r>3. (27) 

Lemma 4 For each a > there exist positive constants K > and /iq 
depending on a and M2 such that, for h < H < ho, hi = h2 = h or 
{/ii,/i2} = {h,H}, and every wl^{-) G Vhi,r and ,wl^{-) £ Vh^^r satisfying 
the threshold condition 

\Ht)-wl,{t)\\.<ahl'\ \\u{t)-wl^{t)\\^<ahl-\ i = 0,l, t£[0,T], 

(28) 

for Wh{t) £ Hl{VL), t £ [0,T], satisfying 

\\wh{t)\\j < 2Qmax(/ii,/i2)^"-', 
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the following bounds hold: 

||-^(^^/i>^^l2)|lo + ll-^(^^i'^'^)|lo - ^W'^hWi, (29) 
\\^(.Wh,wl^)\\_^ + \\j='{wl^,Wh)\\_-^ < K\\wh\\o, (30) 

\\Bh{wh,wl^)\\^ + \\Bh{wl^,Wh)\\(^ < K\\wh\\i, (31) 

||A,^/'(e,(u;^,0)||„ + ||^^^/'(fi,«,u;;,))||Q < K\\wh\\o, (32) 

where T{u, v) can he either (n • V)^; + ^(V • u)v or (n • V)v, and Bh = H/iJ^. 

Proof The proof of the present lemma can be found in that of Lemma 4.4 
in [H] for F{u,v) = {u ■ V)v + i(V • u)v and Wh G Vh.r- With obvious 
changes, the proof is also valid when J-'{u, v) = (n • V)u, as well as when 

Wh Vh,r- □ 

Remark 1 We will apply the above inequalities for Wh = w\_^ — w^^, 
Wh = wj^^ — u and Wh = — u. Let us also remark that the Lemma H] also 
holds if either w]^_^ or is replaced by u. In what follows we will apply 
Lemma [4] to Uh and Vh both satisfying the threshold condition ()28p for an 
appropriate value of a (see |141 Remark 4.1]). 

Lemma 5 For v G {H'^{Q))'^ n V there exists a positive constant K = 
i^(||w||2) such that w € H^iyi)'^ the following hound holds for e = v — w: 

\\A-^U[I'{v,e)+I'{e,v)]\\^^ < K\\v-w\\-i, (33) 

where J-{u, v) can he either (n • \/)v + ^(V • u)v or [u ■ \7)v. 

Proof The proof of this result when J^{u,v) = {u ■ V)v + |(V • u)v can 
be found as part of the proof of [3 Lemma 3.4]. With obvious changes, the 
proof is also valid when J-'{u, v) = {u ■ V)v. □ 

Let us observe that the approximation over the finer grid Uh and the 
recently defined Vh satisfy 

Uh + AhUh + '^hiuH • Vuh) = 'n.hf, Uh{0) = UhUo, (34) 

Vh + AhVh + Uh{u-Vu) = Uhf, Vh{0) = UhUo, (35) 

respectively. Then Ch = Vh — Uh satisfies 

Ch + AhCh + Uh{uH ■ Ve/j) = UhPh,H, ehiO) = 0, (36) 
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where 

Ph,H = UH ■ VVh - U ■ Vu. 

In the proof of Theorem [1] below we wUl use the fohowmg lemmas. 

Lemma 6 Let {u,p) be the solution 0/ ([T|)-([2]) . There exists a positive con- 
stant C such that the following inequality holds for r = 3,4; 

\\A-'U,p,^h\\o < -^^^\\og{H)\H'^\ t G (0,r]. 

Proof Let us write ph^n — Ph h ~^ Ph H-^ where 

Ph,H = {{uH-u)-Vvh), pl^H = {u-y{vh-u)). (37) 
By applying (jl4p we have 

WK^'^hPiuh < Ch^\K^\\, + p-inp{,,||o, j = 1,2. 

To bound ||/oji//||o let us recall Remark [T] and apply to get 

j:^r— 1 

\\Ph,H\\0 < C\\UH - u\\i < C ^(^_.2)/2 ' (-^^^ 

where we have applied (j25p from Lemma [1] in the last inequality. Applying 
([2^ we also get 

\\Ph,H\\0 < C\\vh -U\\i< g ^(^_2)/2 ' (39) 

where in the last inequality we have applied standard bounds for II/j (see 
pT]) ) together with the estimates ([27|) for Zh in Lemma[3l Let us next bound 
||^~^nyojj ji^llo. We will use the decomposition 

Ph,H = ii^ - uh) ■ y){vh -u) + {{u - uh) ■ V)u. (40) 
Then, we obtain 

WA-^Upl^Hh = \\A-^Il{{{uH-u)-V){vh-u))\\o 

+ \\A'^Jl{{{uH-u)-V)u)U (41) 

To bound the second term in ([^T]) we apply ([33l) from Lemma [5] to get 

WA'^n {{{uH - u) ■ V)n) llo < C\\uH - n||_i. 
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Applying then together with (fT2|) and ([27]) we get 

\\UH-U\\-1 < \\UH - Vh\\o + \\vh - u\\-i 

< ^2\'-^mH^^^ + ^H^^\ (42) 

To bound the first term in ([^T]) we argue by duahty, using ([5]), we get 

\\A-^U{{{uH -u)-V){vh -u)) Wo 
= sup (((uH --u) • V)(v/, -'u),A-in(/.) 

110110 = 1 

< sup - u||o||w/i - n||i||^~^n0||oo 

\\<l>\\o=l 

< sup C||u// — ti||o||w/i — "iilli ||^~"^n(/)||2 < C||nH — n||o||fh — 

\W\o=l 

Now, in view of the case r = 2 in ([25]) and using again pT]) and (p7|) we 
conclude 

p-^n (((UH - ^x) • V)K - n)) llo < ^^^H'h^-\ 

Finally, to bound ||y4^^np^ j^||o we apply again (|33|) to bound this norm in 
terms of — n||_i which, as we shown in (IMI), is bounded by Ct(2-r)/2;jr+i^ 
□ 

Lemma 7 Let {u,p) be the solution o/ ([l])-([2]) . There exists a positive con- 
stant C such that the following inequalities hold for r = 3,4; 

\\ph,H\\-i < j^^\log{H)\H^+' + -^^^h^, te{0,T], (43) 

\\Ph,H\\-i < te{0,T]. (44) 

Proof The proof is very similar to the one of the previous lemma. We will 
prove (|13]) since the proof of (jM]) is completely analogous and yet easier. 
We use the decomposition (j37|) . 
For /)| we apply (i30]l to get 

\\pIh\\-i < C\\vh - u\\o < C-r- 



t(r-2)/2 ■ 
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where we have apphed pT]) and (p7|) in the last inequahty. For p\ ^ we will 
use the decomposition ()40p . For the first term in ()40p using ([5|) we have 



= sup {{{uH - u) ■V){vh - u),(j)) 

< sup \\UH - u\\L2d\\Vh- U\\l\\(j)\\^2d/(d-l) 
Il0l|l = l 

< C||ti// — n||i||z;/i — n||i < C 



i(r-2)/2 ^1/2' 



where we have applied ([25]) and (fTTj) and ([27|) in the last inequality. Finally, 
for the second term using ([7]) and (jl2]) we obtain 

\\{{uH - u) ■ V)u\\^i = sup (((li/f - n) • V)n,(/)) 

Il0lli=l 

< sup II^H — 'u||_i||</>Vu||i 



< sup \\UH - u\\-i(\\Vu\\yyi,2d/(d-i)\\4>\\L2d 



+ llVul 



oo 9 I 



f(r-2)/2 I 



□ 



Lemma 8 Let {u,p) be the solution of Then there exists a posi- 

tive constant C such that the discrete velocity Vh defined by (|35p and the 
approximation to u over the finer grid, Uh satisfy the following bound: 

\\A'l\vh{t)-Uhm\o<CH^~^, r>3, / = 0, 1, t e {0,T]. 



Proof Let us consider yh{t) = Af[ eh{t). From (|36]l it follows that 



{t)= I e-^'-'^^'^Aflih{{uH-y)eh)ds+ f e-^'~'^^'^ AfllnPhAs) ds. 
Jo Jo 



Vh 



Applying (|T5|) . and taking into account that as a consequence of (j29j) and 



we have \\A 



(-1+0/2 
h 



I /2 

^h{{uH ■ V)e/i)||o < C\\Ai[ ehWo, it follows that 



\yhit)\\o< n hh^ods + 

Jo yt — s 



e-^'-'^^^J^i^^^p^As)ds 



15 



so that a generalized Gronwall lemma [22, pp. 188-189] allow us to write 



max ||wi,(t)||n < C max 



(45) 



Using (jl6p we obtain 



(46) 



To conclude we apply (I13p and (I44|) from Lemma [7] in the case / = 0, and 
(138)1 and (1391) in the case 1 = 1. □ 



Lemma 9 Let {u,p) be the solution of ([T])-(l2]). Then, there exists a pos- 
itive constant C such that the discrete velocity v^ defined by (jgjp and the 
approximation to u over the finer grid, satisfy the following bound: 

\\A-^^'\vh{t)-Uh{t))h<C\\og{H)\H\ r>3, tG(0,r]. 



Proof The proof follows the steps of the proof of Lemma 4.6 in [14] . Let 
US consider yh{t) = eh{t). From ([36j) it follows that 



yh{t) = I e 

10 



(t-s)AnA-^/2 
Applying (|15p we have that 



A-^^^Uh{{uH-V)eh) ds+ f e-^'-'^^'^A-^/^UhPh,H{s)ds. 

Jo 



' C 
Vt- s 



IVhmo < I -^J\A^''nh{{uH-V)eh)\lds 



+ 



f e'^'-'^^'^A-'^^UhPh,H{s)ds 
Jo 



(47) 



Let us now bound ||yl^ "^n/j((ti// • V)e/i)||Q. We will argue as in the proof of 

[HI (4.23) Lemma 4.4]. Let us first observe that /i^||e/i||i < C|| A^^''^e/i||o. 
Using (fTH) we get 



\A-^UhiiuH-V)eh)\l < Ch''\\{uH-V)eh\\o + \\A-^muH-'^)eh)\\o 

< Ch%h\\i + \\A-^muH-V)eh)\\o 

< C\\A~'/^eh\\o + \\A~'U{{uH ■ V)eh)\\o- 
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Let us now bound the second term on the right hand side above. We write 

\\A-^muH ■ V)eh)\\o < \\A-^m{uH - u) ■ V)eO||o + p-'n((n • V)eOI|o. 

For the first term arguing by duahty we get 

\\A~^U{{{uH - u) ■ V)eh)\\o < C\\UH - u\\o\\eh\\i < CH^eh\\i. 

For the second one, arguing again by duahty and integrating by parts, we 
get 

\\A-^mu ■ V)eh)\\o < C||e;,||_i||^x||2 < C\\yh\\o. 
We finahy obtain 



\\A-'Uh{{uH ■ V)eO||o < C\\yh\\o + CH^ChWi. 
Going back to (jl7|) we obtain 



(48) 



\yh[t)\\o 



< 



c 



\/t - s 



\yh\\ods + 



-{ts)An 4-1/2 



Af^ UhPh,H{s)ds 



+Ct^/^H^ max ||e?,(s)||i, 

0<s<t 





(49) 



so that a generahzed Gronwah lemma ^22^ pp. 188-189] ahow us to write 



Using (|16|) we obtain 



max \\vh{t)\\ri < [ max 
o<t<T"^^ - \o<t<T 



f\-it-s)A,A-y^Yi^p^,jds 
Jo 



+ H max \\eh\\i 



0<t<T 



max^ ||y/i(i)|lo ^ C( ^ I — I niax s 



1 r 

2' 2 / o<sTr" 



^/'^\\Af^'^UhPh,H\L+H'^ max ||e/,||i ), 
" " ' '"J o<t<T J 



(50) 
□ 



where, by applying Lemmas [6] and [8] the proof is finished. 



The proof of the following theorem follows the steps of the proof of |14l 
Theorem 4.7]. 

Theorem 1 Let {u,p) be the solution of There exists a positive 

constant C such that the discrete velocity defined by l\35\i and the approx- 
imation to u over the finer grid, Uh, satisfy the following bound: 



\vh{t) - Uh{t)\\o < :^\log{h)\ {\log{H)\H^) , te{0,T], r>3. 

(51) 
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Proof Let us consider yhit) = t^^'^eh{t). From ([36]) and an easy calculation 

we get 

Vh + AhVh + t^^^UhiuH ■ Vet) = t^/^UhPh,H + ^^^h- 



Then, 

Jo 



+ 



Applying (|32|) we get 

\\A-^^^Uh{uH ■V)eh)\\o < C\\ehh. 
Then, using (|15p we obtain 



(52) 



Jo n Jo 



IVhWo 



ds. 



/o \/r- s 

Applying a generalized Gronwall lemma [22^ pp. 188-189], it follows that 



<C ( max 

\0<t<T 



+ max 

0<s<t 



fe~^^^'~^h"^IlhPh,Hds 
Jo 



-At,,(t-s) fife 
,,1/2 



(53) 



Applying now Lemma [2] and (|16p we have 

max \\yh{t)\L<c{\\og{h)\ max ||s^/2^;;;^n/,p/,,H(s)||n 



0<t<T' 

+ B [ -, - I max IIA. ^^^e/,(s)| 
2 2/ o<t<rll ^ ^1 



where Lemmas [6] and [9] finish the proof. 



□ 



Theorem 2 Let {u,p) be the solution of ([I])-([2]). There exists a positive 
constant C such that the discrete velocity Vfi defined by ()35p and the approx- 
imation to u over the finer grid, Uh, satisfy the following bound for r > 3 



\vH{t) - %(t)||i < j\ \og{h)\ (l \og{H)\H^ + ^ t G (o,r]. 



(54) 
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1/2 

Proof Let us define yh{t) = tA^ eh{t), where eh{t) = Vh{t) — Uh{t). Ar- 
guing exactly as in the proof of Theorem [H instead of (j53|) we now arrive 
at 



max^ ||yh(i)|lo <cf inax^ 



0<t<T 



sAJ UhPhHds 



+ max 

0<t<T 



Applying now Lemma [2] we get 



C\log{h)\( 



max 

\0<t<T 



sAf^^^^UhPh^H 



+ max 

0<t<T 



/I -1/2 



where Lemmas [7] and [9] finish the proof. 



□ 



Lemma 10 Let {u,p) be the solution of ([I])-©- Then there exists a pos- 
itive constant C such that the discrete velocity Vh defined by (fg5l) and the 
approximation to u over the finer grid, satisfy the following hound: 

IIVKW-%W)IIo<c^^', r>4, tG(o,r]. 

Proof Let us consider yh{t) = ^^^e/i(t). From (f36]l it follows that 

yh{t)= [\-^'-'^^'^A-'Uf,{{uH-V)eh)ds+ [\-^'-'^^'^A^'UhPhM^)ds. 
Jo Jo 

We first observe that arguing exactly as in the proof of il4[ Lemma 4.13] we 
get 



+ C^\\eh\\i. 
\/s 



(55) 



Then, 



ivhimo < c 



+ 



■^/|= + \\yh{s)\\o ds 

f e~^'~'^^^Al^IihPh,H{^)ds 
Jo 



Ct^/^H^ max 

0<s<t 



eh{s)\\i 
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Applying now [14^ Lemma 4.9] we get 



max ||wf>(t)||o <C max 
o<t<T"^^ - 0<t<T 



Jo 



+ CH^ max \\eh{t)\\i. 

0<t<T 



(56) 



For the second term on the right-hand-side above we apply Lemma [5J For 
the first one we use the decomposition 

Ph,H = pI,H + Ph,H, Ph,H = ii'^H - u) ■ V{vh - u)) + {{UH - u) ■ Vu), 

We now argue exactly as in [141 (4.60) in Lemma 4.14], replacing one of the 
occurrences of z there by u — uh and making use of (I25p and (I26p with h 
replaced by H. This will allow us to obtain 



max 

0<t<T 



Jo 



which concludes the proof. 



□ 



Theorem 3 Let {u,p) be the solution of ([I])-©. There exists a positive 
constant C such that the discrete velocity Vh defined by (f35j) and the approx- 
imation to u over the finer grid, Uh, satisfy the following bound: 

\\vh{t)-Uhmo<j\log{h)\{\log{H)\H'>), te{0,T], r>4. (57) 



Proof Let us define yh{t) = teh{t), where eh{t) = Vh{t) — Uh{t). Arguing 
as in the proof of Theorem [H instead of (j53p we now arrive at 



max ||y/t(t)||o <C( max 



+ max 

0<t<T 



s^hPh,H ds 



eh ds 



As in the proof of Theorem [21 applying now Lemma [2] to both terms on the 
right-hand side above we get 

^^l^\\yh{t)\\o < C\log{h)\ (^m^ax^||s^-infc/5fc,H||o + Qmax^||^h^e/i|lo) • 



where now Lemmas E] and [TU] finish the proof. 



□ 
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Lemma 11 Let {u,p) be the solution of Then there exists a pos- 

itive constant C such that the discrete velocity defined by (jgjp and the 
approximation to u over the finer grid, Uh satisfy the following bound: 

\\A-'^\vh{t)-Uhm\o<^\logih)\H', r>4, t G {0,T]. 



Proof Setting yh{t) = t^'"^ A^^''^ eh{t) and arguing exactly as in the proof 
of Lemma [H instead of (|49p we now obtain 



< 



c 



\/t- s 



IVhWods + 



-it~s)A.l/2 







+ 



it-s)A,,. -1/2 eh{sl 



s'/'Al'''lihPhMs)ds 



+ CH^ max s^l'^\\eh(s)\\i, 

0<s<t 



SO that a generalized Gronwah lemma [22', pp. 188-189] ahow us to write 



max ||?yh(^)lln <C max 
o<t<T ^ 



+ max 

0<t<T 



0<t<T 
t 

e 



[\-^''^^^-s'/'A-'/'u,p,,H{s)ds 
Jo 







(t-s)A^A-i/2eh{s)_ 



+ CH'^ max t^^'^\\eh{t)\\i 



0<t<T 



Using we obtain 



max \\yh(t]\\f, < CB '. 



) ( max \\sAf^'^UhPh,H\\o + max ||A;^^eh||o) 
/ \0<s<T ' ' II o<s<r J 



+CH^ max t^^'^WehWi- 

0<t<T 



For the first two terms on the right-hand-side above we apply Lemmas [6] 
and 1101 respectively. For the last term we observe that denoting by yh{t) = 
t^l^A]l\h and arguing as in Theorem [T] we get 

max \\yhm\o <c{\ logWI m^x^ \\t^'^ A-^^'^Uhph^Hh 



0<t<T 



0<t<T 



SO that applying now (jl3p and (j44p to bound the first term on the right-hand 
side above, and the case / = in Lemma [8] for the second one, the proof is 
completed. □ 
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Theorem 4 Let {u,p) be the solution of ([I|)-([2]). There exists a positive 
constant C such that the discrete velocity Vfi defined by l\35\i and the approx- 
imation to u over the finer grid, Uh, satisfy the following bound for r > 4; 

\\vH{t) - < ^1 \og{h)\ (l \og{h)\H^ + T'/'^h") , t G (0, T]. 

(58) 

Proof Let yh{t) = t^^'^A^^'^e^ and argue as in the proof of Theorem [2] to 
get 



|yhWllo<C|log(/i)| ( max 



t^/^AV'^'u 



hPh,H 



+ max 

0<i<T 



for t G (0,T]. To bound the first term on the right-hand side above we 
apply ([T3|) and (03]), and for the second one we apply Lemma [TTl □ 



We now summarize the main results of the section in the following theorem. 



Theorem 5 Let {u,p) be the solution of ([I])-©- There exists a positive 
constant C such that the approximation to u over the finer grid, Uh, satisfy 
the following bounds for r = 3,4 and t G (0,T]; 

Mt) - n,(t)||o <^(^l log{h)\\log{H)\H'+' + ^^^h\ 

\\u{t) - u,{t)\\i <^^\log{h)\ (\log{h)\H^+'+T'/'h^)+^^h^-\ 

where in the last inequality we can replace the second |log(/i)| by |log(iJ)| 
in the case r = 3. 

Proof We use the decomposition u — = (u — v^) + {vh — Uh). To 
bound the first term we apply (jlip and (j27p while for the second we apply 
Theorems [HI □ 

Now, we get the error bounds for the pressure. We begin with some 
error estimates for the time derivative of Vh — Uh- 

Lemma 12 Let {u,p) be the solution of ([I])-©. Then there exists a pos- 
itive constant C such that the discrete velocity defined by (jgjp and the 
approximation to u over the finer grid, Uh satisfy the following bound for 
r = 3,4.- 

||t;,(t)-4(t)||-i<^(^|log(/i)|(|log(/i)|/7'-+i + /iO, tG(0,T]. (59) 
In the case r = 3 the second \og{h) can be replaced by log{H). 
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Proof Using ([36|) and taking into account that ||e/i||_i < CH^d^ e/i||o we 
obtain 

||eh||_i < e^llo + ' n,, ((n// • V)eft + llo 

— 1/2 

< llefclli + C||e/j||o + ^hPhM^, 

after using (|5'2p in the last inequahty. Applying now Theorems [T]H] together 
with ([13D and (IMl) we reach (IMl). □ 



The following Lemma is proved in |141 Corollary 4.19] and Proposi- 
tion 3.1 in [23] and P] (see also p] and [E]). 

Lemma 13 Let {u,p) be the solution o/ ([T])-([2|) and let [u^^Ph) o.^'d {vh,gh) 
the approximations defined in 1 1'})- [TS\) and \21\ )- [2^) . respectively. Then, 
the following hound holds for r = 2,3,4 

\\Ph{t)-p{t)\\L2/^<j^^^h'-\ tG(0,T], (60) 
\\9h{t)-p{t)\\L2/^<^^^^h'-~\ tG(0,r]. (61) 

Theorem 6 Let {u,p) he the solution of There exists a positive 

constant C such that the approximation to p over the finer grid, ph, satisfies 
the following hound for t S (0, T] and r = 3, A: 

\mt) - p{t)\\Lyw < ^jSm^""' + ^(^5)72 1 (I iog(MI^''+' + h^) ■ 

In the case r = 3 the second log(/i) can he replaced by log{H). 
Proof We use the decomposition 

\\Ph -pWl^/r < \\Ph - 9h\\L^m + \\9h -p\\ 



To bound the second term on the right-hand-side above we apply ()6ip . For 
the first one subtracting (f2T]l from (fTUI) and applying the inf-sup condition 
pOj) we obtain 

P\\Ph - ghh^/R < \\uh - Vhh + \\{uH ■ V)eh\\-i + \\ph,H\\~i + \\eh\\~i- 

We first observe that applying (j30|) we get \\{uh • V)eh\\-i < C'||e/i||o. To 
bound ||/9/i^//||_i we apply ()i3|) . Now, the proof concludes applying Theo- 
rems [THU together with (|59p . □ 
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We state in the following theorem the results that can be obtained for the 
mini-element. 



Theorem 7 Let {u,p) be the solution of ([I|)-(I2|)- There exists a positive 
constant C such that the approximations over the finer grid computed using 
the mini-element, {uh,Ph), satisfy the following bounds for t S (0,T]; 

||n(t) -u(t)||o < CH^ + Ch\ 
\\u{t)-u{t)\\i < C\log{h)\H^ + Ch, 
\\Phit)-pit)\\L2/R < C\logih)\H^ + Ch. 

4 Semi-discretization in space. The second two- 
grid algorithm. 

As in the previous section we will concentrate on the approximations ob- 
tained with the Hood- Taylor mixed finite element and r = 3 or r = 4 and 
we will state at the end of the section the results that can be obtained for 
the mini-element method with a much simpler analysis. 

In the second algorithm we consider, the first level, as before, is given 
by the standard mixed finite-element approximation to ([I])-©, that is, the 
solution of (fT7|) - (fT8]) with initial condition uniO) = ^niuo)- In the second 
level we solve a linearized problem on a finer grid and given n/i(0) = HhUo, 
we compute Uh{t) ^ Xh,r and Phit) G Qh,r-ii t £ (0,r], satisfying, for all 
0h G Xh,r and Tph G Qh,r-i 

(uhAh) + (yuh,V(l)h) + b{uH,Uh, (ph) + b{uh,UH, 4>h) + {Vph, (ph) = 

{f,4>h) + b{uH,UH,(l)h), (62) 

(V-n/„Vh) = 0. (63) 

Observe that the approximation uh is the result of one step of Newton's 
method for the Galerkin {uh,Ph) approximation in Xh^r x Qh,r-i (equa- 
tions (fT7|) ~ (fT8]) with H replaced by h) with {uh,Ph) as initial approxima- 
tion. For this reason, in this section we study the error eh = Uh — Uh- 
It is easy to obtain that 

eh + Aheh + BhiuH,eh) + Bhiet, uh) = ^hPh,H, (64) 

where 

Ph,H = —F{eh,H,eh,H), 
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where, here and in the sequel, 



The analysis in this section is closely related to that in the previous 
section. However, some extra results are needed. We shall use the following 
two bounds, 

||</'h|li2./(.-i) < c UhWl'^ \\A]l^4>h\\]l\ y4>h G Vh,r, (65) 
Uhho^ <C\\Al^^cl)h\\l^^\\Ah(l)h\\o\ y(l>h(^Vh,r, (66) 

which follow from Corollary 4.4 and Lemma 4.4 in [23]. Also we shall use 
the following two bounds 

\\^h'^^^Bh{eh^H,eh,H)\\o < C\\eh,H\\l^^\\(^h,H\f/'^ , (67) 

(68) 

with C independent of h, and where here and in the sequel Bh{vh,Wh) = 
IlhF{uh,Wh)- Both are easily obtained by duality arguments, the first one 
from (j65p and the second one from (|66p . Notice also that as a conse- 
quence of (ISSI-dMl) and ([IS]) we have that lle-*^''^/,!!^^ < Ct~^/^ \\(j)h\\o 

and ||e~*^''A^/^(^/j||^2d/(d-i) < Ct~^^^ \\4'h\\o so that by using duality argu- 
ments together with these two inequalities the following two bounds easily 
follow 

\\e-('--^^'^Bh{eh,H,eh,H)t <j^^^^\\eh,H\\o\hM\r (69) 
||e-(*-^)^''-4^'i?h(6.,i^,6,,^)||,<^^-^||6,,H||^ (70) 

Lemma 14 There exists a positive constant C = C{M2) such that 

\\u-uh\\loo < CH^'^. 



Proof We will use the fact that, due to Lemma 4.3 and 4.4 in (23j . and 
CoroUary 4.4 in [23]. 

\\VIiHu\\L.<C\\Au\\, (71) 

We write u — uh = {I — IIh)u + {Hhu — uh)- Applying Q, we have 

~ '^hWq ^ CH^^"^, where in the last in- 
equality we have applied pT]) and ([25]) . On the other hand, applying 
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(4.43)] 



11(1 - UhM^^ <C\\{I- Uh)u\\]!,^ II V(/ - UH)ufl' 

<C\\iI- UhMY^ (l|Vn||^6 + ||Vn;,n||^6 



a/2 



Now, where, in the last inequahty we have apphed ([5]) and |23t Lemma 4.4]. 
Furthermore, applying pT]) . ^ and (j7T]) the proof is finished. □ 

Lemma 15 Let {u,p) be the solution of Then there exists a positive 

constant C such that the approximations Uh and Uh to the velocity u over 
the fine mesh satisfy the following bound: 

\\A',l\uh{t)-Uh(t))\\o<CH'/^-', r>3, / = 0,1, t e {0,T]. (72) 

Proof Follow the arguments in the proof of Lemma[8]to obtain (j46p . Now, 
for / = 1 we write s^/^ \\UhPh,H\\o < Cs^l'^ ||eh,H|lioo II eh,// Hp so that apply- 
ing Lemma [T3] and (|25p the proof of the case / = 1 is finished. For / = 0, 
applying (f671) we have 



.1/2 



S (-ys ||e/i,i/||o 



l|3/2 



so that applying 



11 11 i|\l/2||l/2 II 

= U|e/i,i/|lo l|e/i,//|lij IP' e/i,i/||]^, 
the proof is finished. □ 



Lemma 16 Let {u,p) be the solution o/ ([T])-([2]) . Then there exists a positive 
constant C such that the approximations Uh and Uh to the velocity u over 
the fine mesh satisfy the following bound: 

\\A-^\uh{t)-Uh{t))h<CH\ r>3, tG(0,r]. 

Proof Follow the arguments in the proof of Lemma [TOl but notice that 
now due to the terms (e^ • V)uh and (V • UH)eH, instead of ([55]) we have 

e-(*-)A.^-in,((^^ . v)e,)|[ <C (-^ + -L) p-ie,||o 



+ ^ {—^\\eh\\i + ^lle/illo 



Thus, instead of (156]) we arrive at 



max ||yh(t)||o <C' max 



0<t<T 



f e'^'''^^'^A^'Uhph,H{s)ds 
Jo 

+ C(H^ max \\eh(t)\\i+H'^ max \\eh 

^ 0<t<T 0<t<T 



0<t<T 
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Thanks to Lemma [15] we have that the last two terms on the right-hand side 
above are bounded by CH^^/'^. For the first one, applying first (j68p and 
then (|26p we conclude that it is also bounded by CH^ . □ 



Lemma 17 Let {u,p) be the solution o/ ([T])-([2]) . Then there exists a positive 
constant C such that the approximations and Uh to the velocity u over 
the fine mesh satisfy the following bound: 



Af^'/\uh{t) - Uhimio < -^H\ r>3, t e {0,T]. 



C 



(73) 



Proof Let ynit) = t^^'^^h ^^"^^hit) and follow the arguments in the proof of 



Lemma [9] so that instead of (j47p we now have 

rt c 



Wh 







< 







s^l^{Bh{eh,UH) + Bh{uH, 



+ 



+ 



^/t^s 

\-^'~'^^>^Al^'^Iins^/^pn,H{-^)ds 



ds 







(74) 



Now observe that by using ||n/i-(s) — "^(5)11^ < CH^ ^ / s^^'^ ^ instead of 
we now have 

\Al\s^/\BH{eh,UH) + BiXuH,eh)\^<c(^\yh{s)\\Q+H^\\eh\\^+H^\\^^^^^ 
Thus, instead of ([5U|) we now get 

<CB ( mj^x^4A-in;,p;,,H||o + max^||A-ie/,(s)| 



max % 



+ C{H^ max lle/i 111 + i/^ max ||e/i||o). 

^ Q<t<T 0<t<T ' 





(75) 



Due to Lemma [15] we have that the last two terms on the right-hand side 
of ()75p are o(H^), and due to Lemma [T6l the second one is 0{H^). Finally 
due to (f68l) the first one can be bounded by C||s^/^e/i^i^ ||Q||s^/^e/i^//||-^, which, 
due to ([25D is also 0{H^). ' ' □ 

Theorem 8 Let {u,p) be the solution of ([I])-©. Then there exists a pos- 
itive constant C such that the approximations Uh and Uh to the velocity u 
over the fine mesh satisfy the following bound: 

C 



\uh{t)-Uh{t)\\Q<-H\ tG(0,r], r>3. 



(76) 
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Proof Let yhit) = tehif) and argue as in the proof of Theorem [T] so that 
similarly to (j53p we now get 



max % t) n 



<c 



, max 

\o<t<T 



+ max 

0<t<T 



Jo 



-^i^^^-'^Ch ds 



Using (j69p to bound the first term on the right-hand side above, and ([16 
for the second one, we get 



max \\yhmo<c(T^/^ max ||fi/2e;,,^(t)|| ||tV2 

0<t<T ^ ''lu y n^+^'pii ' 'MUM > ^ 'mi 

+ B 



0<t<T' 

-,-] max||ti/2A-^/V(t)| 



so that applying (j25p and Lemma [T71 the proof is finished. 



□ 



Theorem 9 Let {u,p) be the solution of ([I])-©. Then there exists a pos- 
itive constant C such that the approximations and Uh to the velocity u 
over the fine mesh satisfy the following bound: 



C 



||n,(t)-n,(t)||i<^(^-^F'-+\ te{0,T], r = 3,4. (77) 

Proof Let yhit) = t^'-'^^'^ A]!"^ eh{t) and follow the arguments in the proof 
of Lemma [8] so that now, instead of ()45p we get 



max 

0<t<T 



\yK{t)\l <cf max f e-^'-')^>^ s^'-^^'^ A]I^IihPh.H{s) ds 

\0<t<T Jq 



+ 



(r-1) 



(78) 



Applying (f70|) to bound the first term on the right-hand side above and (fTBI) 
for the second one, we have 



max 

0<t<T 



|y.(t)|lo<c(rV4 maX^||il/2,,_^(i)||J|t(-2)/2^^^^(^)||^ 



+ CB 



0<t<T 

-,-] max||t(^-2)/2e;,(t)|l 



(79) 



Due to (j25p the first term on the right-hand side above is bounded by 
CH^H^^^ = CH^'^^. For r = 4, the second one is bounded in Theorem [8l 
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When r = 3, we may write = pe/j(t)||Q''^ l|e/i(t)||o^^ so that ap- 

plying Theorem [5] and Lemma [T^ the second term on the right-hand side 
of ^ is bounded by CH^/^H'^/^ = o{H^) □ 

Remark 2 For r = 3 it is possible to prove the bound 

\\uh{t) - un{t)h < J {h^^^ \log{h)\ + H'''/') , t e (0,r], r = 3. 

To do so, apply Lemma[2]and (j67p to bound the first term on the right-hand 
side of (|78p and the same bound as before for the second term. 

Finally, repeating (with obvious changes) the analysis in Section [3] for 
the pressure, the following result is easily proved 

Theorem 10 Let {u,p) be the solution of There exists a positive 

constant C such that the approximation to p over the finer grid, ph, satisfy 
the following bound for t G (0, T] and r = 3, 4; 

\\ph{t)-Ph{t)\\LyR < j^SmH^"-', t G (o,r]. (80) 

We now summarize the main results of the section in the following theorem. 

Theorem 11 Let {u,p) be the solution of ([l])-([2]). There exists a positive 
constant C such that the approximations {uh,Ph) satisfy the following bounds 
for r = 3,4 and t e (0,T].- 

\Ht)-u,{t)\\o<jH' + j^h^. 

Proof We use the decomposition u — u = {u — Uh) + {u^ — u) and apply 
(|25p to bound the first term and Theorems [8] and [9] for the second. For the 
pressure, using the decomposition p — ph = (p — Ph) + (Ph —Ph) and applying 
()60p and Theorem [10] the proof is finished. □ 

Finally, with a much simpler analysis, that we do not detail here for brevity, 
the following result can be proved 
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Theorem 12 Let {u,p) be the solution of ([II)-©. There exists a positive 
constant C such that the approximations over the finer grid computed using 
the mini-element, {uh,Ph), satisfy the following bounds for t G (0,T]; 

\\uh{t) - u{t)\\o < CH^ + Ch\ 
\\uh{t) - u{t)\\i < CH^ + Ch, 
\\Ph{t)-pmL2/M < CH^ + Ch. 

5 Fully discrete case. 

In this section we consider the fuhy discrete case. Let us assume that we 
integrate in time equations ()17til8p and (jl9H20p for the first method or equa- 
tions (jlTtilSp and (162II63P for the second method using the backward Euler 
method or the two-step backward differentiation formula (BDF). In the case 
of the two-step BDF method the first step is carried out using the backward 
Euler method. We will denote by (U^,P^) the fully discrete Galerkin ap- 
proximations to the velocity and pressure at the time level tn = nk for 
< n < AT and /c = At = T/N. We wih denote by (C/^,^^) the fully 
discrete approximations to the velocity and pressure over the finer grid at 
the time level tn- 

Let us denote by = unitn) — and by e^ = Uh{tn) — UJ^ the 
temporal errors in the approximations and UJ^ respectively. Let us 
denote by tt^ = pnitn) — Ph ^"^^ by vr^ = Ph{tn) — PJ^ the temporal errors 
in the approximations and P^ respectively. In [15] we have proved the 
following error bounds. There exist constants Ci^ and ko such that for k < ko 
the temporal errors of the Galerkin approximation satisfy the following error 
bounds 

l|eS-||o + inP/fe^||o < Qol^' lkHllL2(n)/R < Qo (2;o-i)/2 ' 1 < < ^> 

tn tn 

where = 1 for the Euler method and Iq = 2 for the two-step BDF. Let 
us remark that using ||e^||i < C||A]|^^e^||o < C||e^|lo^^||^_H'e^|lo''^ error 
bounds in the norm are also obtained in a straightforward manner. 

Using the same technique developed in [H] it can also be proved that 
analogous error bounds hold for the approximations over the finer grid. More 
precisely, for both the first and second algorithms there exist constants Cig 
and ko such that for k < ko the temporal errors of the two-grid approxima- 
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tion satisfy the following error bounds 

\\eh\\o + tn\\Aheh\\o < Cioj^pr^ \\^h\\L'2{n)/R < Qo (2;o-i)/2 ' 1^"^-^' 

*n tn 

where = 1 for the Euler method and Iq = 2 for the two-step BDF. 
Finally, using the decompositions 

<in)-Uh = Htn) - Uhitn)) + el, 
p{tn)-Pi: = {p{tn)-Ph{tn)) + TTl 

the error bounds of the fully discrete approximations are obtained as the 

sum of the spatial errors (the errors in the semi-discrete approximations we 
have already bounded in the previous sections) plus the temporal errors. 
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